Scaling properties of the asymmetric exclusion process with long-range hopping 
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The exclusion process in which particles may jump any distance I > 1 with the probability that 
decays as / - ( 1+<T ) is studied from coarse-grained equation for density profile in the limit when the 
lattice spacing goes to zero. For 1 < a < 2, the usual diffusion term of this equation is replaced 
by the fractional one, which affects dynamical-scaling properties of the late-time approach to the 
stationary state. When applied to an open system with totally asymmetric hopping, this approach 
gives two results: first, it accounts for the er-dependent exponent that characterizes the algebraic 
decay of density profile in the maximum-current phase for 1 < a < 2, and second, it shows that in 
this region of a the exponent is of the mean-field type. 

PACS numbers: 02.50.Ey, 05.40.-a, 05.60.-k 



I. INTRODUCTION 

The asymmetric exclusion process (ASEP) is one of the 
simplest models for transport in which the net flow of par- 
ticles is maintained by contact with two reservoirs at dif- 
ferent densities ph and pr. It emerges in the wide range 
of non-equilibrium phenomena, like biological transport 
0, 0, surface growth 0, 0] and traffic flow In the 
presence of the net current flowing through the system, 
the boundary conditions generally play an important role 
in determining the bulk properties [6|. The phenomeno- 
logical domain-wall approach Q and the exact solution 
|l [§] revealed that the model exhibits phase transitions 
both of the first and the second order. Description of 
these new phenomena contributed to considerable efforts 
taken in the last several decades in order to determine 
the nature of phase transitions in system held far from 
equilibrium. Important aspect of these efforts concerns 
a better understanding of those parameters (e.g. sym- 
metry of the Hamiltonian, dimensionality and range of 
the interactions) that usually determine an underlying 
universality among various models of equilibrium phase 
transitions, but remain important for out-of-equilibrium 
phase transitions as well (for a recent review on univer- 
sality classes in non-equilibrium lattice models, see |10l|). 

In that spirit, one of striking features of ASEP is 
the robustness of its phase diagram to various modifi- 
cations. Among a number of different extensions of this 
model that have been proposed, the universal character 
[EEL 03 of the continuous phase transition has been ver- 
ified in a number of cases, including ASEP with parallel 
update [HI, [HI, partially asymmetric exclusion process 
[la ] and particle- wise disorder With the purpose 

to examine the impact of the long-range hopping on the 
character of the phase transitions, a generalized model 
was introduced [17] in which particles may jump any dis- 
tance I > 1 with the probability that decays as l~( 1+a \ 
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Due to the long range of hopping, the exchange of parti- 
cles with reservoirs is possible at each site of the lattice, 
which can be compared to the system with a bulk reser- 
voir [lfl [l^| . For a > 1 , this generalized model has the 
same phase diagram as the short-range case consisting of 
the low-density, high-density and the maximum-current 
phase, but different effects are found at the transition 
lines. Besides the localization of the domain- wall at first- 
order phase transition, the continuous phase transition to 
the maximum-current phase is accompanied by the ex- 
ponent that differs from the short-range value 1/2 in the 
region 1 < a < 2, where its dependence on a was given 
by the conjecture based on numerical simulations [17] . 

In the present work we show that the conjectured ex- 
pression for the er-dependent exponent can be obtained 
in the mean-field approximation. We consider the hydro- 
dynamic approach that gives coarse-grained equation for 
density profile and apply it to an open system. Within 
the same formalism we examine dynamical scaling prop- 
erties of this model in more general asymmetric and sym- 
metric cases with periodic boundary conditions and con- 
firm the predictions by numerical simulations. 

The work is organized as follows. In Section II we 
consider the process on the infinite lattice and find the 
continuous limit of lattice equations in a general case 
including both the symmetric and the asymmetric ex- 
clusion process. In Section III we analyze the effects of 
introducing the boundary conditions that correspond to 
those of the maximum-current phase and deduce the an- 
alytical expression for the cr-dependent exponent. Direct 
numerical solution of lattice equations in the mean-field 
approximation is also given. The Section IV is dedicated 
to the analysis of the relaxation to the stationary state. 
A brief summary of results is given in Section V. 



II. HYDRODYNAMIC APPROACH 

Let us first consider the exclusion process on the infi- 
nite one-dimensional lattice, where each site n is either 
occupied by a particle (r„ = 1) or empty (r„ = 0). Dy- 
namics of this process is described as follows. At any 
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given time t, a randomly chosen particle at site n at- 
tempts to jump either to the site n — I with a probability 
q, or to the site n + I with a probability p = 1 — q. Dis- 
tance I > is chosen according to the probability distri- 
bution pi = ;~( 1+<T ) /£(<r + 1), where C,(z) is the Riemann 
zeta function. Generally, one has p ^ q (the asymmet- 
ric exclusion process) , or p = q (the symmetric exclusion 
process). In the limit a — > oo the hopping reduces only 
to nearest neighbors, and the process is identical to the 
exclusion process with short-range hopping. 

This dynamics corresponds to a continuous-time 
Markov process in which the probability P(C, t) of the 
system being in a state C = {r n } ne z at time t evolves 
according to the following master equation 



C)P(C',t) - J2 W(C -> C')P(C,t), 
c 

(1) 

where W(C — > C) is the transition probability per unit 
time for a system to go from the state C to the state C . 
In the model considered here, the non-diagonal elements 
W(C — C) are equal to q • p n _ m if C = C<- n ^ and 
m < n, or to p ■ p m -n if C = C^ m,n ^ and m > n, where 
(j(n,m) j g ^- ne configuration obtained from C by moving 
a particle from a site n to an empty site m (if possible) . 
Starting from the master equation ([T]), lattice equation 
for the time-dependent average density profile (T„)(t) = 
J2 C T n P(C,t) reads 



<r n )(f) = <X«), 



(2) 



where 



"(P - <l)^2Pr [(1 - Tn)A+T„ + T n A~T n ](3) 



r>0 



and the following notation has been used, A+r„ = r n+r — 

T n and A r T n = T n T n — r . 

Basically, hydrodynamic equation that we are inter- 
ested in can be obtained by a suitable coarse-grained pro- 
cedure when the lattice spacing a goes to zero. This is 
usually referred as the hydrodynamic limit and describes 
the time evolution of a system at large time and space 
scales where the stochastic details of the particular pro- 
cess have been smoothened out. So far, rigorous results 
have been obtained in the short-range [2ll. |22| and re- 
cently in the long-range case for p = q [27| . In the short- 
range case, the equation for the coarse-grained density 
p(x,t) is given by the inviscid Burgers' equation in the 
asymmetric case p^ij, 



and by the normal diffusion equation in the symmetric 
case p = q, 



dp 
dt 



i crp 



(•5) 



In the long-range case |27|, diffusion equation (J5J) is re- 
placed by space-fractional diffusion equation in the region 
1 < a < 2, 



— = ^A a (j)(x,t), l<cr<2, 



(6) 



where A^ is fractional Laplacian (|A.5[) and v a = 
-2 P r(-cr)cos(7ra/2)/((a + 1) > 0. 

To study the asymmetric case we adopt the non- 
rigorous approach that consists of taking the mean-field 
approximation (r n r TO ) — » (T n )(r m ), n ^ m and applying 
it directly to lattice equations (T2J) and Notice that 
the same approach, when applied to the short-range case, 
gives the correct hydrodynamic limit although the start- 
ing mean-field assumption is only approximate (in the 
symmetric case, the assumption is not necessary since 
equation is linear). The reason for this lies in the fact 
that for certain initial profiles stationary state becomes 
a product measure (i.e. a factorized state without corre- 
lations), so in this case the mean-field approximation be- 
comes exact. It is worth mentioning that the factorized 
state is preserved even upon the introduction of long- 
range hopping *28], provided that a hopping probability 
has a finite mean (in case pi ~ Z - ( 1+CT ', this is true for 
a > 1). 

Applying the mean-field approximation to lattice equa- 
tions ([5]) and ©, we obtain the following equations 



= £f(A r + ^-A^ n ) + 

r>0 

+{Ap + (f> n )( P - g)£p r (A+0„ + A-0„)(7) 

where <f> n (t) = (T n )(t) — p denotes deviation from the uni- 
form profile of density p and Ap = p — 1/2 is introduced 
in order to distinguish two cases, p — 1/2 (Ap — 0) and 
p ^ 1/2 (Ap 7^ 0). To obtain equation for the macro- 
scopic profile 4>(x, t) in the continuous limit, we adopt 
the procedure given in (24[. Basically, one assumes that 
4> n (t) are coefficients in the Fourier series of some func- 
tion 4>(k, t) defined on [-K/2, K/2] 



(8) 



where x n = na and a — 2w/K is lattice constant. Then, 
( 4 ) Eq. © can be written in the Fourier space 
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dt 



4>{k,t) = (f)(k,t)[D(ka)-D(a)]+Ap4>(k,t)B(ka) + 

x rK/2 i-K/2 

+ 1F2 / dM(fci>*)0(fc2,t) 

K J-K/2 J-K/2 



■ e l[kl+k2 ~ k)na B{ka) 

n— — oo 

where D(ka) and B(ka) are given by 



D{ka) = - [Li a+1 {e lKa ) + Li a+l {e 



-ika \ 



(9) 



B(ka) = (p - g)[L V+ i(0 - Lv + i(e- ,fe ")] (10b) 



and Li v {z) is the polylogarithm function. Using the se- 
ries representation of the polylogarithm function valid for 
(10a) non-integer v ^ 1, 2, 3 . . . [29|, one obtains 

I 



D(ka) -D(0) = 



1 



2C(a + l) 



2r(-a)c 03 -|fcra- + 2 y c(<t + \t 

V ; 2 1 1 ^ (2n) 



^(zfc) 2 "a 2 " 



(11a) 



B(ka) 



p-q 

C(^ + l) 



-2ir(-a)sin^-sgn(k)\k\' T a a + 2 ^ ^L. 

n=l 



2nl 



(2n-l)! 



2n-l 2n-l 



r 



(lib) 



Finally, one defines the macroscopic density profile 4>(x, t) 
as the inverse Fourier transform of <f>(k, t), obtained from 
cj)(k,t) in the limit a — > after the appropriate scaling 
t — > f/a z has been taken, where z is the lowest exponent 
in a in (|1 la|l and (|llb[) . The latter should not be confused 
with dynamical exponent z considered in Section IV. 

The symmetric case (p — q). From (jTTaj) it follows 
that z = min{a,2}. For a > 2 under diffusive scaling 
[z = 2) one obtains the normal diffusion equation 



dt 



d 2 (b 

"2-^-2, cr > 2, 
OX 



(12) 



with the diffusion coefficient v-i = £(er — l)/2£(<7 + l) > 0. 
On the other hand, in the region 1 < a < 2 one obtains 
the space-fractional diffusion equation 



1< a < 2, 



(13) 



where v a = — T(— a)cos(TTa/2)/((a + 1) > 0, as in Ref. 

da. 

The asymmetric case (p ^ q). From (jllap and (jllbp 
it follows that z = min{a, 1}. For <j > 1 under Eulerian 
scaling (z = 1) this gives the inviscid Burgers' equation 
with the additional drift term —vd/dx<j>(x,t) 



dt 



dx 



K(j) 



dx ' 



a > 1, 



(14) 



where the collective velocity v and k are given by 



v = {p-q)(l-2p)X(a) 



dj( P ) 



dp 



p=p 



K=-2(p-q)X{a) 



d 2 j( P ) 



dp 2 



(15) 



(16) 



p=p 



In the above expressions, A(<r) is the average hopping 
length and j(p) is the macroscopic current 



j(p) = (p-q)X(a)p(l-p), X(a) = 



CM 



C(<7 + 1) 



(17) 



The original inviscid Burgers' equation (01 is then recov- 
ered either by Galilean transformation x — > x — uf or by 
taking p = 1/2. 

We can go further and take a look at the next higher- 
order terms in (|lla[) and (|llb[) . For 1 < a < 2 there are 
two terms, both proportional to a" . The first one follows 
from (|llap and corresponds to the fractional Laplacian 
A a (j) (also known as the Riesz fractional derivative), while 
the other one is nonlinear and reads (f>H a (j), where H a is 
defined in (|A.7|) . For a > 2, the next term in piap 
corresponds to the usual diffusion term term A</>(x) and 
is proportional to a. If we neglect the nonlinear term 
4>Ha(j), we obtain following equations 
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at 



= av 2 A(f) a 

at 



dx 



1< a < 2, (18) 



ox 



<j > 2, 



(19) 



where <fi(x,t) is replaced with cj) a (x,t) to emphasize its 
dependence on a. These additional viscous terms do not 
represent the true lowest-order correction, but are rather 
kep t due to the smoothening (regularizing) effect (see 
[3l| and references therein) they impose on the solutions 
of the inviscid Burgers' equation that could otherwise 
develop discontinuities (shocks) [23|. As we show in Sec. 
IV, they also play an important role in selecting the way 
in which the system relaxes to the stationary state. 



III. APPLICATION TO AN OPEN SYSTEM 



We now wish to apply hydrodynamic equations (|18[) 
and (fTTH) to a finite system in contact with left and right 
reservoirs of densities pl = a and p^ = 1 — (3, re- 
spectively. Furthermore, we consider only the station- 
ary limit where d(p(x,t)/dt — so that <j>(x,t) — > (f>(x). 
This case was already studied by Monte Carlo simula- 
tions with random sequential update [l7| in the totally 
asymmetric case (p = 1, q = 0) for various values of 
a and (3. The obtained density profiles reproduced the 
phase diagram of the short-range model provided that the 
current was renormalized by the average hopping length 
\ L {a) = a(<?)/U(<J + 1), where ( L (z) = £f =1 r<l+>) is 
partial sum of the Ricmann zeta function £(z) and L is 
the number of sites on the lattice. 

However, a difference was observed at the first-order 
transition line and in the maximum-current phase. Apart 
from the localization of a domain-wall at the first-order 
transition line a = (3 < 1/2, the exponent that char- 
acterizes the algebraic decay of density profile in the 
maximum-current phase appeared to be a-dependent. In 
particular, the results obtained by Monte Carlo simula- 
tions in the maximum-current phase for system of size 
L showed that a deviation of density from its bulk value 
1/2, Ap(n,L) = |(r n ) — 1/2|, obeys the scaling relation 



Ap(n,L)=L-»f(n/L) 



(20) 



where f(x) 
be 



for x <C 1/2 and \x was conjectured to 



P 



a-1 1 



(21) 



In the rest of this Section, we show that this particular 
dependence on a can be explained within the hydrody- 
namic approach. Our approach follows the one of Krug 



|{|, who applied the viscous Burgers' equation (fT9f in 
the stationary limit to a finite system of size I in con- 
tact with two reservoirs in order to investigate boundary- 
induced phase transitions in the asymmetric exclusion 
process with short-range hopping. Imposing the bound- 
ary conditions p(0) — po and p(l) = 0, he was then able 
to show the emergence of characteristic length £, diverg- 
ing for po greater than the critical density p* = 1/2 and 
finite £ ~ (p* — po) -1 for po < p* . Thus the continuous 
transition occurs between the phase in which p{x) decays 
exponentially to the bulk value p ^ p* and the power- 
law in which £ diverges and p{x) displays algebraic decay 
with the exponent p, = 1, p(x) — p* ~ a/x. This ap- 
proach gives the correct qualitative picture that explains 
the second-order phase transition in the asymmetric ex- 
clusion process with open boundaries, but the exponent 
itself is wrong and should be 1/2 @, Q instead of 1 [32| . 
The reason lies in the fact that in general the steady state 
of an open system is no longer given by the product mea- 
sure but displays correlations, while these arc by default 
neglected in the mean-field approximation. 

Let us start by assuming that particles jump only in 
one direction, so that p — 1 and q = 0. Due to the 
finiteness of the system, a distance 1 < I < L is now 
chosen according to the probability pi = l~ t - 1+(T > /Cl{c + 
1). As it was pointed out in our previous work [17j . 
the long-range hopping in a finite open system raises the 
necessity of "non-local boundary conditions" , in the sense 
that the exchange of particles with reservoirs now takes 
place at each site as if there was a bulk reservoir [l8l Il9l . 
[20| . As the Fig. 1 shows, one can see that, for example, 
a particle at site n is removed from the lattice with the 
probability [3 n , which is the sum of probabilities of all the 
possible jumps outside the lattice. This includes all the 
jumps of size I for which L — n + I < I < L, multiplied 
by the probability 1 — pn, = (3 that the right reservoir is 
empty 



Pn 



P 



E 



1) t *+f 

1 l = L-n+l 



(22) 




PR = 1-P 



FIG. 1: (Color online) A schematic picture of an open system 
with L = 6 sites. A particle at site n may jump to the right 
reservoir in a number of ways for which n + I > L (dashed 
lines), all adding up to the total probability j3 n to remove the 
particle from a system. Similarly, a particle may be added 
from the left reservoir to an empty site n only from those 
sites in the left reservoir that are within a distance n < I < L 
from the site n. 

Similarly, a particle from the left reservoir is added to an 
empty site 1 < n < L with a probability a n , which is the 
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total probability of all the possible jumps from the left 
reservoir to an empty site n, multiplied by the probability 
Pl = a that the left reservoir is not empty 



L 



cr+l ' 



(23) 



A precise meaning of these boundary conditions can be 
given by noticing that in case of the open boundary 
conditions, the lattice equations for the average density 
(t„)(<) in the mean-field approximation 



where the most dominant part of the expression (|28|) 
follows from the "boundary" condition <fi(x) — <fi(0) for 
x < 0. [As the nonlinear term (pHv^i is concerned, one 
can show that it is of the order 0(a cr_1+ ' t x~ a ~f i )] Since 
the left hand side of Eq. (p~8j) must be zero in the sta- 
tionary regime, terms (f2"5)) and (f2"T|) must be of the same 
order, which gives the scaling exponent p = (cr — l)/2, 
both in powers of a and x. 

For a > 2, the boundary effect due to the fixed value 
of 4>(x) for x < requires that one includes both aA(j)(x) 
and a a ~ x A a (j>(x) in the equation for 4>(x). This gives the 
following estimate for aA(f>(x) 



|(n) = a x (l - ( n » - J2 Pm -i(n)(l (r m » - /h(n), 

(25) 



n-l 



— (r n ) = a n (l - (r n )) + ^ p n _ m (r m )(l - (r n )) 



XI Pm-n<r„)(l - (7Yn» - /3„(r n }, (24) 



m— n+l 



m— 2 



a^ 2 A</)(a;) ~ 1 1 A* ( 1 + n)a 



(29) 



remain the same. 
min{a — 1,1 + /i} that 



while the expressions (|27|l and 
This leads to equation 2/2 
has two different solutions depending on the value of cr: 
p = (cr - l)/2 for 2 < a < 3 and p = 1 for a > 3. To 
summarize, we can write 



/x = mm 



cr - 1 



(30) 



=a L (l-(r L ))+> PL-m{T m )O—{TL))-0L(TL), 



at ^— ' 

rn— 1 

(26) 

can be written in the same form as in {2} and ([3]) if we 
extend the system for additional L sites to the left and 
the right, but require that (r„) = a if — L < n < and 
= 1 — /3ifi<n< 2L. [A similar reasoning that 
fixes the values of function in the extended region was 
used for the numerical solution of the boundary-value 
problem that involves fractional derivatives [33j .] 

In order to examine the scaling property (|20p , we 
suppose that away from the left boundary (but sub- 
stantially far away from the right one), density profile 
4> n = \{r n ) — 1/2| decays algebraically with some un- 
known exponent p, > 0, <f> n ~ nT^ = cfi jx^ . A rough 
estimate of the length scale beyond which this asymp- 
totic behavior sets in is given by Is ~ a/(po — l/2) 1 /^, 
where po = (n). To employ the hydrodynamic equation 
for this problem, we assume that the size of the system is 
large enough that the influence of the right reservoir may 
be ignored. In that case, we may fix the value of <p{x) to 
be equal to 4>(0) = po — 1/2 for all x < 0. Inserting 4>{x) 
in the Eq. (fT5)) yields the following estimate of various 
terms for x > Is 



K0[X) 



( d(j){x) 
dx 



-PIk^x- 2 *- 1 , 



(27) 



l v a A a (j)(x) 



aC,{a + 1 



cr— 1 cr 

-a X 



-0(a 



For 1 < cr < 2, this result is exactly the same as (|21j) . 
However, it fails to give the correct value of a for which 
the short-range regime sets in: according to (|30"f this 
value is a — 3, instead of a = 2. This failure arises as 
a result of neglecting the correlations in the mean-field 
approach, so that for a > 2 the non-local effect of the 
boundaries is overestimated. 



Numerical solution of discrete mean-field 
equations 



(28) 



The expression (|30| can be checked directly by the nu- 
merical solution of the stationary lattice equations 
(|26|) in the mean-field approximation. This case reduces 
to a problem of finding a zero of a system of L nonlinear 
equations in L variables, which can be done numerically. 
For this purpose, we used the HYBRD algorithm taken 
from the MINPACK library 34] for various a, (3 and a. 
In the region 1 < a < 2, the results reproduce the phase 
diagram well and the profiles coincide with the results of 
the Monte Carlo simulations (Fig. [2] and [3]). The only 
exception is the line a = j3 < 1/2 with a sharp domain- 
wall located in the middle (Fig. |4j. This is similar to 
the mean-field solution in the short-range case [H| which 
does not take into account fluctuations of the position of 
the domain wall. These fluctuations for the long-range 
case [l7| can be taken into account using the domain- 
wall approach in a same manner as it was done for the 
short-range case with the bulk reservoir [35| . where the 
domain-wall performs random-walk in a potential well. 
As far as the maximum-current phase is concerned, den- 
sity profiles satisfy the scaling relation (f2"0")) along with the 
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mean- field exponent p, for all 1 < a < 2 (Fig. [5]), but the 
profiles themselves no longer match those obtained in the 
Monte Carlo simulations as soon as a > 2 (Fig. [6|. This 
is expected since for a > 2, p, reads min{(a — l)/2, 1} 
while the correct value is 1/2. 
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FIG. 2: (Color online) A comparison of density profiles ob- 
tained in the mean-field approximation (line) and by Monte 
Carlo simulations (symbols) for a — 0.2, j3 — 0.7 and a — 1.8 
(low-density phase). 
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FIG. 3: (Color online) A comparison of density profiles ob- 
tained in the mean-field approximation (line) and by Monte 
Carlo simulations (symbols) for a — 1.0, ft = 1.0 and a — 1.8 
(maximum-current phase) . 

Let us conclude this Section by giving a simple pic- 
ture that accounts for the cr-dependent exponent fi for 
1 < a < 2 and its change to the short-range value 1/2 
for a > 2. First, notice that the nonlinear term <pd<p/dx 
appears in the hydrodynamic equations no matter what 
the range of the hopping is. This implies the possibil- 
ity of reproducing the same tr-dependent exponent by 
replacing the long-range hopping with the short-range 
one, provided that the nearest-neighbors hopping rate 
is increased by the factor X(a) and that the bulk reser- 
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FIG. 4: (Color online) A comparison of density profiles ob- 
tained in the mean-field approximation (line) and by Monte 
Carlo simulations (symbols) for a = f3 = 0.2 and a = 1.5 
(coexistence line). 
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FIG. 5: (Color online) Deviation <j> n (L) of a density profile 
from its bulk value p = 1/2, obtained from the numerical 
solution of mean-field equations for various system sizes L — 
800, 1600, and 3200 (a = = 1.0) and for a = 1.2 and 1.8 
(from top to bottom). The profiles (f>„(L) for the same a are 
scaled to the profile cj>„(Lo) {La = 3200) according to Eq. 
(|20[) with exponent /i = fi given by (1301) . 



voir is kept with the same a-dependent rates a n and (3 n . 
Indeed, density profiles obtained by Monte Carlo simu- 
lations of such a modified model reproduce the scaling 
property (|20|) with the exponent (a — l)/2 for 1 < a < 2 
and 1/2 for o > 2. As far as the large-scale behavior 
is concerned, the original model is thus reduced to the 
short-range one with the additional external field that 
originates from the long-range exchange of particles be- 
tween the reservoirs and the bulk, while all the other 
contributions that are result of the long range of hop- 
ping are of higher order. [For example, this is the case 
with the hopping from one site in the bulk to another. 
As already mentioned earlier in this Section in context of 
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FIG. 6: (Color online) A comparison of density profiles ob- 
tained in the mean-field approximation (line) and by Monte 
Carlo simulations (symbols) for a — 1.0, /3 = 1.0 and a — 2.5 
(maximum-current phase) . 



hydro dynamic approach, in the maximum-current phase 



the nonlocal term 4>H a <p is of order 0(a° 



-l+fi~,—a—fi 



) and 



does not affect the scaling properties.] Since the exter- 
nal field itself exhibits a power-law dependence in the 
position on the lattice, it affects the scaling exponent in 
the maximum-current phase. However, the latter is true 
only for 1 < a < 2, where the influence of the boundaries 
is stronger than that of the correlations. The boundary 
value <T = 2 may be reasoned by looking at the average 
distance 7l(ct) = J2 n n ' <W a = J2 n ( L ~n+l) ■ /?„,//? 
at which the particles are created and annihilated [17| . 
Since 7l(c) diverges with the system size L as L 2 ~ a for 
1 < a < 2 and tends to a finite value for a > 2, the 
influence of the external field becomes localized near the 
boundaries for <j > 2 and the exponent /j, changes to its 
short-range value 1/2. 

At this point, it is useful to give a more detailed com- 
parison between the present model and the (short-range) 
model of TASEP with Langmuir kinetics [ij, Hi]. In the 
latter case, particles are created and annihilated in the 
bulk with rates £Ia and f2o, respectively, that depend 
on the system size as L^ 1 . In the more general case 
[20| . these rates have been extended to Qa,^d ~ L~ a , 
1 < a < 2, which is the same power- law dependence 
displayed by a n and (3 n [ItJ with a = a. As a result 
of this dependence, both models exhibit localization of 
the domain- wall at the coexistence line a = /3 < 1/2 
for 1 < a < 2, where the width of the domain- wall 
£-a/2 j g determined only by the a-dependent potential 
well in which the domain-wall performs a random walk 
[l7l I20L [3o] | . On the other hand, the distinction between 
these two models becomes pronounced in the maximum- 
current phase, where the fact that the rates a n and (3 n de- 
pend on position on the lattice becomes important. This 
is in contrast to the short-range case with bulk reser- 
voir where the creation and annihilation of particles with 
homogeneous rates does not affect the scaling exponent 



which remains 1/2 for all 1 < a < 2, but merely deter- 
mines the characteristic time that particles spend on the 
lattice (H. 



IV. DYNAMICAL SCALING 

Another issue we wish to address here concerns the 
dynamical-scaling properties of the late-time approach 
to the stationary state. In the short-range exclusion pro- 
cess, the longest relaxation time r scales with a system 
size L as r ~ L z where z equals 2 for the symmetric 
(p = q ) (36| and 3/2 for the asymmetric case (p =/= q) 
[36l . [13, l38l. |39|| . These exponents were brought into con- 
nection with dynamical exponents of the corresponding 
hydrodynamic equations with additional noise term. In 
particular, dynamical exponents z — 2 and z — 3/2 cor- 
respond to those of the noisy Edwards- Wilkinson (EW) 
[4C| and the noisy Burgers' equation [41(, respectively, 
where the latter can be mapped to the Kardar-Parisi- 
Zhang (KPZ) equation [26| of surface growth. 

Generally, the late-time characteristics of a system can 
be probed by looking at the two-point autocorrelation 
function which in a translationally invariant system takes 
the form 



C(:M) = (0(O,O)#M)>, 



(31) 



where the averaging (...) is taken over the noise histories. 
For processes described by EW or KPZ equations, C(x, t) 
is known to give the following scaling relation, 



C(x,t) = x 2x ~ 2 F(t/x z ), 



(32) 



where \ denotes the roughening exponent and is equal 
to 1/2 in both cases (ill. |42||. 

Both EW and KPZ equations have been generalized 
by replacing the usual diffusion term with the fractional 
one [23, 53| • This gives the fractional EW equation 



dh A , 



< a < 2, 



(33) 



where <j)(x,t) = dh(x,t)/dx and r){x,t) is the noise. A 
simple dimensional analysis giyes the values of exponents 
X = (c — l)/2 and z = a [43|]. On the other hand, the 
fractional KPZ equation, 



dh A (dh 

m=» A ° h+ 2(dx 



< a < 2, 



(34) 



displays a more complex behavior with various values of 
X and z depending on the value of a and whether the spa- 
tial correlations of the noise are relevant or not [25] . Par- 
ticularly, in case these correlations are irrelevant, there is 
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a weak-coupling regime for a < 3/2 with the same criti- 
cal exponents as in the fractional EW equation (f3"3"|) , i.e. 
X = (it — l)/2 and z — a, and a strong-coupling regime 
for a > 3/2 with the usual KPZ exponents \ — 1/2 and 
z = 3/2. 
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FIG. 7: (Color online) Time-decay of the two-point autocor- 
relation function C(0, t) ~ t~ 1//z for the symmetric exclusion 
process with long-range hopping on the half-filled periodic 
lattice (L = 10 4 , averaged over 10 7 independent MC runs). 
Dashed lines refer to the expected value of z = min{a, 2}. 




t 

FIG. 8: (Color online) Time-decay of the two-point autocorre- 
lation function C(0, t) ~ t _1 ^ z for the asymmetric exclusion 
process with long-range hopping on the half-filled periodic 
lattice (L = 10 4 , averaged over 10 7 independent MC runs). 
Dashed lines refer to the expected value of z = min{a, 3/2}. 

To compare this picture with the late-time approach to 
the stationary state of the exclusion process with long- 
range hopping, we computed the corresponding two-point 
autocorrelation function C(i — j, t) using Monte Carlo 
simulations in a system with periodic boundary condi- 
tions. The averaging was taken over 10 independent 
runs with L = 10 4 for times up to t — 100 Monte Carlo 
steps/site. Dynamical exponent z was then extracted 



from the time-decay of the peak C(0,t) ~ t~ x l z p3 |. 
The results for various a are presented in Fig. [7] and 
Fig. [5] for the symmetric and the asymmetric case, re- 
spectively. In the symmetric case, dynamical exponent is 
in a good agreement with z — min{a, 2} predicted by the 
fractional (1 < a < 2) and the normal (a > 2) EW equa- 
tion. In the asymmetric case, results presented in Fig. [H] 
confirm a change in the exponent at a = 3/2 from the 
weak-coupling to the strong-coupling regime according 
to the predicted value z = min{a, 3/2} of the fractional 
KPZ equation with irrelevant spatial correlations of the 
noise. In the interpretation of the Ref. (25[, this change 
reflects a tendency of a system to relax through its fastest 
"component" that has a lower value of z. 



V. CONCLUSION 

In this work we studied the large-scale dynamical and 
stationary properties of the exclusion process with long- 
range hopping where each particle may jump any dis- 
tance / > 1 with the probability that decays with I as 
Pi ~ l~( 1+cr '. We started on an infinite lattice with no 
boundary conditions present and calculated the contin- 
uous limit of lattice equations in the Fourier space [24| . 
For the symmetric hopping this yields the same result 
as in Ref. 27} in which the fractional diffusion equation 
replaces the usual one for 1 < a < 2. For the asymmet- 
ric hopping, we adopted the non-rigorous approach that 
consists of decoupling lattice equations by means of the 
mean-field approximation. As in the short-range case, 
the result is still given by the inviscid Burgers' equa- 
tion, but the particle current j(p) = (p — q)X{a)p{\ — p) 
has increased by the factor X(a) — £(°')/C( cr + 1). A 
true signature of the long range of hopping was found 
by inspecting the lowest-order corrections in powers of a 
which in this case correspond to the fractional Laplacian 
for 1 < a < 2 and the usual one for a > 2. If these terms 
are kept, the equation becomes equivalent to the deter- 
ministic part of the fractional (1 < a < 2) [25[ and the 
original KPZ equation (cr > 2) [2^]. To check this connec- 
tion we analyzed the late-time approach to the stationary 
state in a finite system with periodic boundary condi- 
tions. For that purpose we extracted dynamical expo- 
nent from the time-decay of the autocorrelation function 
obtained from Monte Carlo simulations and compared it 
to the dynamical exponent z = min{a, 3/2} of the frac- 
tional KPZ equation [25[ . A very good agreement of these 
two exponents lead us to the conclusion that the system 
always relaxes through its fastest "component" that has 
the lower value of z [25[ . Similar analysis was also carried 
out for the symmetric case yielding a very good agree- 
ment with the known exponents z = a and z — 2 of the 
fractional (1 < a < 2) and the usual Edwards- Wilkinson 
equation (a > 2), respectively. 

In case of the open boundary conditions, the above 
approach was applied only phenomenologically. In ad- 
dition, the long range of hopping in a finite system in- 



9 



troduces "non-local boundary conditions" by means of 
the inhomogeneous external field that creates and anni- 
hilates particles at each site. In the maximum-current 
phase, this causes the change in the exponent fi that de- 
termines the algebraic decay of density profile. The re- 
sulting exponent obtained in the mean-field approxima- 
tion justifies the earlier conjecture [l7| that \i = (a— l)/2 
for 1 < a < 2. We argue that the short-range regime sets 
in immediately after, for a = 2. 



where f(k) — T{f(x)}. If we now take the linear combi- 
nation of -ocT^x and x T^%o given by 



A CT /(x) = - 



2cos(tt(t/2) 



(A.5) 



we get the Riesz fractional derivative (also known as frac- 
tional Laplacian) with the following property 



APPENDIX: RIESZ FRACTIONAL DERIVATIVE 



f{A a f(x)} = -\k\°f(k). 



(A.6) 



Consider the left- and the right-hand Riemann- 
Liouville fractional derivatives a T>^ and X T>% of order a, 
respectively, whose action on the suitable function f(x) 
is defined as 130] 



aKf(x) 



1 



r(n — er) dx n J a 



n — a—1 



(A.l) 



On the other hand, the first term in (jllb[) corresponds 
to the following linear combination of -ooT^x an( i x^Zo 



H a f(x) = 



_ -oa Ly x x^o 



2sin(na/2) 



(A.7) 



with the following property 



'*™=wh<&L f ^- x) > (A - 2) 

where n is the smallest integer exceeding a. For a = — 00 
and b = 00, integrals (jA.lj) and (|A.2j) have the following 
simple property with respect to the Fourier transform 



f{H <T f(x)} = -isgn(k)\k\°f(k). 



(A. 



Both A CT and H a are special cases of more general 
fractional derivative defined as the inverse of the Feller 
potential (sec [30] for details). 



F{xVZof(x)} = (ik) u f(k), 



(A.3) 



(A.4) 
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